Search for: All records

We present the discovery of PSR J1947–1120, a new huntsman millisecond pulsar with a red giant companion star in a 10.3 day orbit. This pulsar was found via optical, X-ray, and radio follow-up of the previously unassociated gamma-ray source 4FGL J1947.6–1121. PSR J1947–1120 is the second confirmed pulsar in the huntsman class and establishes this as a bona fide subclass of millisecond pulsars. We use MESA models to show that huntsman pulsars can be naturally explained as neutron star binaries whose secondaries are currently in the “red bump” region of the red giant branch, temporarily underfilling their Roche lobes and hence halting mass transfer. Huntsman pulsars offer a new view of the formation of typical millisecond pulsars, allowing novel constraints on the efficiency of mass transfer and recycling at an intermediate stage in the process. 

Abstract PINTis a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework withinPINTto characterize the single-pulsar noise processes present in pulsar timing data sets. This framework enables parameter estimation for both uncorrelated and correlated noise processes, as well as model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated data sets as well as a real data set of PSR B1855+09. We also describe the new features implemented inPINTsince it was first described in the literature. 

Abstract We present the discovery of a new optical/X-ray source likely associated with the Fermi γ -ray source 4FGL J1408.6–2917. Its high-amplitude periodic optical variability, large spectroscopic radial-velocity semiamplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of 0.027 ± 0.003 M ⊙ . Considered jointly with the more massive redback millisecond pulsar binaries, we find that the “spider” companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses of between ∼0.07 and 0.1 M ⊙ . X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient “redback-like” shocks. Together black widows and redbacks dominate the census of the fastest spinning field millisecond pulsars in binaries with known companion types, making up ≳80% of systems with P spin < 2 ms. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical 1.4 M ⊙ , and many of the highest-mass neutron stars claimed to date are black widows with M NS ≳ 2.1 M ⊙ . Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses. 

Abstract Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 ± 0.8 mm s^–1yr^–1, which deviates from the acceleration predicted by equilibrium Milky Way (MW) models at a 4σlevel. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the MW interacting with orbiting dwarf galaxies (∼1 mm s^–1yr^–1), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in Gaia DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 60 kyr, a companion mass of 0.3M_⊙(indicative of a white dwarf or low-mass star), and a semimajor axis of 1900 au. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. 

Abstract We test the impact of an evolving supermassive black hole mass scaling relation (M_BH–M_bulge) on the predictions for the gravitational-wave background (GWB). The observed GWB amplitude is 2–3 times higher than predicted by astrophysically informed models, which suggests the need to revise the assumptions in those models. We compare a semi-analytic model’s ability to reproduce the observed GWB spectrum with a static versus evolving-amplitudeM_BH–M_bulgerelation. We additionally consider the influence of the choice of galaxy stellar mass function (GSMF) on the modeled GWB spectra. Our models are able to reproduce the GWB amplitude with either a large number density of massive galaxies or a positively evolvingM_BH–M_bulgeamplitude (i.e., theM_BH/M_bulgeratio was higher in the past). If we assume that theM_BH–M_bulgeamplitude does not evolve, our models require a GSMF that implies an undetected population of massive galaxies (M_⋆≥ 10¹¹M_⊙atz> 1). When theM_BH–M_bulgeamplitude is allowed to evolve, we can model the GWB spectrum with all fiducial values and anM_BH–M_bulgeamplitude that evolves asα(z) =α₀(1 +z)^1.04±0.5. 

Abstract Neutron Star Interior Composition Explorer has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. Two model parameters track different types of background events, while the third is used to predict a low-energy excess tied to observations conducted in sunlight. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4–12 keV; 50 detectors) of 0.87 c s −1 , but in 5% (1%) of cases, the rate exceeds 10 (300) c s −1 . Model residuals persist at 20%–30% of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c s −1 (3 σ , single GTI) at 0.4–12 keV, equivalent to 3.6 × 10 −12 erg cm −2 s −1 for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c s −1 at 0.3–2.0 keV, or 4.3 × 10 −13 erg cm −2 s −1 for a 100 eV blackbody. These limits would be four times lower if exploratory GTIs accumulate 10 ks of data after filtering at the level prescribed for faint sources. Such filtering selects background GTIs 85% of the time. An application of the model to a 1 s timescale makes it possible to distinguish source flares from possible surges in the background. 

Abstract Pulsar timing array observations have found evidence for an isotropic gravitational-wave background with the Hellings–Downs angular correlations between pulsar pairs. This interpretation hinges on the measured shape of the angular correlations, which is predominantly quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipolesℓ≥ 2, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and nonzero at the ∼95% confidence level and a Bayes factor of 200. When we include multipolesℓ≤ 1, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz, which has also been noted in previous analyses of the NANOGrav 15 yr data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations. 

Abstract Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio frequency. However, pulse arrival times are also affected by radio-frequency-dependent “chromatic” noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (NANOGrav) and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar data set. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from ${\log }_{10}{A}_{\mathrm{RN}}=-{14.1}_{-0.1}^{+0.1}$to $-{14.7}_{-0.5}^{+0.3}$, and the spectral index broadens from ${\gamma }_{\mathrm{RN}}={2.6}_{-0.4}^{+0.5}$to ${\gamma }_{\mathrm{RN}}={3.5}_{-0.9}^{+1.2}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. 

Abstract Magnetars, isolated neutron stars with magnetic-field strengths typically ≳10 14 G, exhibit distinctive months-long outburst epochs during which strong evolution of soft X-ray pulse profiles, along with nonthermal magnetospheric emission components, is often observed. Using near-daily NICER observations of the magnetar SGR 1830-0645 during the first 37 days of a recent outburst decay, a pulse peak migration in phase is clearly observed, transforming the pulse shape from an initially triple-peaked to a single-peaked profile. Such peak merging has not been seen before for a magnetar. Our high-resolution phase-resolved spectroscopic analysis reveals no significant evolution of temperature despite the complex initial pulse shape, yet the inferred surface hot spots shrink during peak migration and outburst decay. We suggest two possible origins for this evolution. For internal heating of the surface, tectonic motion of the crust may be its underlying cause. The inferred speed of this crustal motion is ≲100 m day −1 , constraining the density of the driving region to ρ ∼ 10 10 g cm −3 , at a depth of ∼200 m. Alternatively, the hot spots could be heated by particle bombardment from a twisted magnetosphere possessing flux tubes or ropes, somewhat resembling solar coronal loops, that untwist and dissipate on the 30–40 day timescale. The peak migration may then be due to a combination of field-line footpoint motion (necessarily driven by crustal motion) and evolving surface radiation beaming. This novel data set paints a vivid picture of the dynamics associated with magnetar outbursts, yet it also highlights the need for a more generic theoretical picture where magnetosphere and crust are considered in tandem. 

Abstract The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency,ν) nor a change in dispersion measure alone (which would produce a delay proportional toν⁻²). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015.